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Abstract 

Continuous-time random walks are generalisations of random walks frequently 
used to account for the consistent observations that many molecules in living cells 
undergo anomalous diffusion, i.e. subdiffusion. Here, we describe the subdiffusive 
continuous-time random walk using age-structured partial differential equations with 
age renewal upon each walker jump, where the age of a walker is the time elapsed 
since its last jump. In the spatially-homogeneous (zero-dimensional) case, we follow 
the evolution in time of the age distribution. An approach inspired by relative 
entropy techniques allows us to obtain quantitative explicit rates for the convergence 
of the age distribution to a self-similar prohle, which corresponds to convergence to a 
stationnary profile for the rescaled variables. An important difficulty arises from the 
fact that the equation in self-similar variables is not autonomous and we do not have 
a specihc analyitcal solution. Therefore, in order to quantify the latter convergence, 
we estimate attraction to a time-dependent “pseudo-equilibrium”, which in turn 
converges to the stationnary profile. 
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1 Introduction 


1.1 Brief model description 

Recent methodological advances in cell biology allowed the measurements of the displacement of 
single molecules (or assemblies thereof) in living single cells. Those investigations have consis¬ 
tently reported that the random displacement of molecules inside cells often deviates from Brow¬ 
nian motion, with the mean squared displacement that does not scale linearly with time, 

as in Brownian motion, but sublinearly, with a power-law behavior: oc [71 El [T5115]. 

This behavior is usually referred to as “anomalous” diffusion or “subdiffusion”, since /i < 1 
usually, for non-active transport (for a review see e.g. [H]). 

Continuous-time random walks (CTRW) are one of the main mechanisms that are recurrently 
evoked to explain the emergence of subdiffusion in cells. CTRW were introduced fifty years ago 
by Montroll and Weiss as a generalisation of random walks where the residence time (the 
time between two consecutive jumps) is a random variable r with probability distribution (fiir) 
(see m for a review). If the expectation of r is defined, for instance when r is dirac-distributed 
or decays exponentially fast, one recovers the “normal” Brownian motion. However, when the 
expectation of r diverges, for instance when (()(r) is heavy-tailed, ())(r) oc with 0 < /r < 1, 

the CTRW describes a subdiffusive behavior, with (r^(t)) oc 

One great achievement of CTRW is that they can readily be used to derive mean-field equa¬ 
tions for the spatio-temporal dynamics of the random walkers. Indeed, starting from (/>(t), 
combinations of Laplace and Fourier transforms lead to a “subdiffusion” equation for the den¬ 
sity of random walkers located at position x at time t: dtp{x,t) = D where 
is a generalised diffusion coefficient and 'DI~^ is the Riemann-Liouville fractional deriva¬ 
tive operator jHIII]. Such a fractional dynamics formulation is very attractive for modelling 
in biology, in particular because of its apparent similarity with the classical diffusion equation. 
However, contrarily to the diffusion equation, the Rieman-Liouville operator is non-Markovian. 
This non-Markovian property becomes a serious obstacle when one wants to couple subdiffusion 
with chemical reaction IHIIIHIE]. 

Here, we take an alternative approach to CTRW that maintains the Markovian property of 
the transport equation at the price of a supplementary independent variable. We associate each 
random walker with an age a, that is reset when the random walker jumps. In one dimension 
of space, we note n(t, cc, a) the density probability distribution of walkers at time t that have 
been residing at location x during the last span of time a. The dynamics of the CTRW is then 
described with an age-renewal equation with spatial jumps that reinitialise the age: 

{ dtn(t, X, a) + dan(t, X, a) +/3(a)n(t, X, a) = 0 , t > 0, a>0, xSK 

n(t, a;, a = 0) = /3(a')uj(x — x')n(t, x', a')dxda' (1) 

n{t = 0, X, a) = n°(x, a). 

The kernel to describes the spatial distribution of jump destinations (typically a Gaussian 
distribution centred at the origin position), and the function /3(a) gives the jump rate. Since we 
are mostly interested here in the subdiffusive case (where the expectation of the residence time 
diverges), we will focus throughout this article on the case: 

a/3(a)- )-/rG(0,1). (2) 

a—>oo 

The precise meaning of the limit will be given later on. The limit p in eq. 0 is the subdiffusion 
exponent: for p, > 1, eq.Q describes a diffusive process, whereas for 0 < /r < 1 the mean 
time a particle has to wait between two consecutive renewals diverges and the mean squared 
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displacement exhibits subdiffusion with exponent /x. The distribution of residence time (fiir) 
evoked above is related to the jump rate as: (/)(r) = /3(r) exp (— /J" /3(s) ds). Note that this age- 
structured approach is not uncommon in the CTRW literature [mil]. Our main contribution 
here is to use it in conjunction with approaches borrowed from the study of partial differential 
equations. 

In the present article, we restrict our attention to the temporal evolution of the age distribu¬ 
tion of the walkers. To this end, we simplify the problem by considering its spatially-homogenous 
version, namely: 


{ dtn{t,a) + dan{t,a) + I3{a)n{t,a) = 0 , t > 0, a>0 

n{t,a = 0)= P{a')n{t,a')da' (3) 

n{t = 0, a) = nP{a). 


1.2 Self-similar solutions 

The only steady state solution of eq.([^ in is 0, which doesn’t allow us to describe the dynamics 
of the system in a satisfactory way. Hence the search for self-similar solutions. An educated guess 
is that they should be of the following form, with A(t) to be determined: 

n(i,a) = -j|^w( ln(l^-b t)^ ,a/A{t)y 

T b 

Let us consider, for the sake of simplicity, an initial condition supported on [0,1). By injecting 
the previous expression into eq. (§, we hnd that the natural choice A(t) = 1 + t preserves the 
initial condition n(0, 6): = n°(&) = w{0, 6): = w°{b), and yields : 

( drW + db{{l — b)w) + e'^fi{e^b)w = 0 

< w{t, 0) = e'^/3{e^b)w{T,b)db ( 4 ) 

[ w ( 0 , b) = w°{b). 

Note that for an initial condition supported on [0, A+), A{t) = A+ -|-1 is a better choice and 
leads to a similar analysis. If the initial condition is not compactly supported, the tail of the age 
distribution can influence the convergence rate we give below. 

It is important to note that the previous system is not autonomous, for the term e'^P{e'^b) 
depends on b. This rescaling does not lead to a classical steady state, and we could not find 
a particular solution of the previous equation. However, we may look for a stationary state 
satisfying formally the following equation, since we consider here /3(a) ~ 

5h((l - 6)IToo) + y VLoo = 0 

b 


where the boundary condition cannot be stated as an equality since Woo is expected to blow up 
at 0, but can be understood as an equivalence as e tends to 0 of Woo(£) and /“ yIToo(fe)d6. 
This leads us to define the self-similar equilibrium as: 


Woo{b) 


CcxD 

bf^{l - 


(5) 


which is called the arcsine distribution, or Dynkin-Lamperti distribution. Coo is defined such that 
||kLoo||i = 1. Under some conditions, we can expect that w{T,b) will converge to eq.([^ when 

T —> 00. 
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A similar result in probability theory appears in Feller’s book [5] tome II, chapter XI, espe¬ 
cially in section 5 and onwards, where the renewal problem is tackled by considering the waiting 
time before the n**' renewal. For an introduction to renewal theory, see the eponymous chapter 
(8.6) in [1]. However, no convergence rate is given for our infinite mean waiting time problem in 
any of these books, and we have been unable to locate such a convergence rate in the subsequent 
literature. Recent developments in Ergodic Theory for mildly related problems (see chapter 
8.11 of [T] for an introduction to Darling-Kac theory), have yielded convergence rates, that are 
optimal in certain cases, as shown in |10j and m- 


1.3 Main results 

Throughout the article, the following set of hypotheses will intervene. Hypothesis (HI) will be 
used in properties of convergence without a rate while hypothesis (H2) will allow convergence 
rate estimates. 


(HI). P is a positive, bounded, and non-increasing function satisfying 


lim a/3(a) = ^ € (0,1). 

a— 


Remark. We will always assume P to be non-increasing for the sake of simplicity (in particular 
in theorem]^ propositionsand 20 and lemmas 13 and 22). The monotonicity can be replaced 
by the following hypothesis: 

/3 is a positive, bounded function satisfying lima_>.oo aP{a) = /i € (0,1), such that P defined 
as follows 

/3(a) = inf /3(a:) (6) 

£cG[0,a] 


also satisfies lima_>.oo ap{a) = yt. 

This leads to minor changes in the proofs, the loss of a multiplicative constant in the affected 
results and replacing P hy P where it corresponds. 


(H2). P satisfies (HI). Additionally, P{a) = 
such that 


u 

\-\-CL 


l5(s)|ds < 


-\-g{a), where g G and there exist K,a > 0 
K 

(1-f a)“' 


Remark. For the sake of clarity we will investigate separately the particular case 5 = 0, called 
the ’’reference case”. Then, all our results will be extended to the general case at the expense of 
the convergence rates. 

Due to the specific shape of W^o and to the boundary condition, it is difficult to investigate 
in a direct way the evolution of ||w — IFoo||i: the methods we describe subsequently fail to 
do so. However, we could recover a quantitative explicit convergence rate with respect to a 
“pseudo-equilibrium” W which will be proved to converge in to IFoo. 


Definition 1. We define the pseudo-equilibrium W over IR+ x [0,1) as follows : 


W (r, b) 


C(r) 

eB(e-b)(l _ 6)1-/^' 


(7) 


where B(a) = P(s)ds and C is defined so that ||IF(t, OIIli = 1- 
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In particular in the reference case Pia) = it may be written as 


W{t, b) 


C{T)e->^^ 

(e-^ + b)t^{l - by-^^' 


Note the similarity between this expression and that for II^oo in eq.([^. In the following, we 
obtain explicit convergence rates of w(t, b) to W, developing proofs based on Relative Entropy 
estimates. Importantly, we show that W converges to W^o at the same rate (up to multiplication 
by a constant) as the rate with which w converges to W. Hence, the convergence to the pseudo¬ 
equilibrium W yields a very good estimate of the convergence to the self-similar equilibrium 
Woo- Finally, we carry out Monte-Carlo simulations of zero-dimensional CTRW to illustrate and 
question the optimality of our main analytical results. 


Definition 2. For the moment and for the sake of simplicity, let us define: 


■H(r) = ||u;(t,-) - lE(r, •)||ii([op]). 

Remark. Later on, we shall define more generally % as & relative entropy, the distance being 
a particular case more suited to our purposes. 

In this paper, we prove the following propositions: 

Proposition 3. Under hypothesis (HI), we have: 




Our first quantitative result is a convergence rate for the reference case of hypothesis (H2). 
Theorem 4. Let P{a) = Then we have the following convergence rates: 

W)5^1(14=.)“] +<=-“-'‘>'[3^1 V 

^ (H(0) + 8r) 



H-e 


A modified, yet analogous, convergence rate still holds for g O'. 

Theorem 5. Suppose hypothesis (H2) holds. 

If a > 1 — fi, we recover the optimal rate of convergence 

fA:(e -^^if fj.^ 1/2 


nr) < 


Kre 


= 1 / 2 - 


If Oi < 1 — p,, we need to distinguish between several cases: 


K{e 


-ar I 




nr) < < 


if — fj, 

K{l + T)e~>^^, if a = fj. < 1 — fj., 

-ke-^^), if a = I- 1/2, 

iL(l + r2)e-^/2^ 


if a = p, = 1 — p, = 1/2. 


if ti=\. 

( 8 ) 
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We finally reinterpret our results in terms of non-rescaled variables, for instance in the 
reference case /3(a) = 

Corollary 6. Assume nP is supported in [0,1) and j3{a) = then if we denote 

iVoo(t,a) = a<l+t, 

jo, a > 1 + t. 


Then if p-^ 1/2, there exists K such that 


K 


\n(t,.) A^oo(t, .)||i < + 


K 


If p = 1/2, then we have 


|n(f,.) -iVoo(t, .)lli < 


{1 + ty {1 + ty-f^' 

K{l + log{l+t)) 


\/l + t 


Remark. An analogous version in non-rescaled variables can be given for theorem 


1.4 Outline of the paper 


The paper is organized as follows. In Sectionwe set the entropic structure of the equation and 
the main properties of the pseudo equilibrium W. In particular we establish (non quantitatively) 
that 


lim 

r—>-oo 


||u;(t, •)-W^(t,-)||i = 0, 


lim 

r—>-oo 


||W^oo(r,-)-H"(r,-)lli =0, 


proving thereby 


lim \\w{t, •) 

T—^OO 


W^oo(t,-)||i =0. 


Section deals with quantitative convergence rates towards the pseudo-equilibrium W, proving 
Theorems 1^ and A convergence rate for Htu — lToo||i, some effects of initial conditions on the 
convergence rates, and convergence rates in non-rescaled variables are dealt with in Section [^ 
Finally, we show the results of some simulations in the eponymous section. 


2 Entropic structure 

Even if we are mainly estimating norms, we see our proof as a specific case of relative entropy 
inequalities. Rates could be obtained following the lines of our proofs for other entropies. 

2.1 contraction for compactly supported solutions 

The first evidence of an attractor is the contraction of compactly supported solutions. If 
we consider two initial data supported in [0,1), non-negative and of mass 1 and the 

associated solutions wi,W 2 , then we have the following property (we take /3(a) = for this 
computation) 

|wi(t, b) — W2 {t, b)\db < D{t) 
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where 


D{t) = 


b) - W2(t, b))db 


e-^ + b 


\wi{t, b) -W2 {t, b)\db. 


Since mass is conserved i.e., wi — W 2 = 0, we have easily 


D(t) = 




M 


e”’" + b e-'^ + 1 
1 


(wi(t, b) - W2iT, b))db 


< 




[ -z:l^^\wi{T,b) - W2{T,b)\db 
Jq e +0 

\wi{t, b) - W2 {t, b)\db 




e”"^ + b e~'^ + 1 
/■i 




|wi(t, b) -W2 {t, b)\db. 
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Thereby, we obtain 


^ lo ~ 'W2{T,b)\db < - J \wi{T,b) - W2{T,b)\db. 

And this leads to ^ 

[ \wi{T,b)-W2{T,b)\db = 0{e~>^^). 

Jo 

In the next section we identify the attractor towards which solutions converge. 

2.2 Pseudo equilibrium 

We start by recalling the definition of what we call the pseudo equilibrium. We recall Definition 

m 


W (t, b) = 


Cir) 


(1 - &)l-MeS(e-b) 

where B{a) = ld{s)ds and C is defined so that ||1 D(t,- lUii = 1. 

Remark. By definition, 

IT(r,& = 0) = C'(t). (9) 

Firstly we establish the fact that W{t^ ■) is an approximation of W^o 
Lemma 7. Assume hypothesis (HI). Then, defining the Dynkin-Lamperti distribution as in 

(§.■ 

- bY-^ 

Wooib) = 7T 


/g b~^^{l — b)^^~^db 


have 


lim ||IF(r, •) - bFoolli = 0. 

r=+<50 


Proof. We start with the model case f3{a) = In this case, we can write 

wi-.t)- (o-'+o-'-o-i'r-' 

+ b)-t^(l - b)f^-^db 

and the result is immediate. 

For the general case, we use the following useful bound on W: 
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Lemma 8. Under hypothesis (HI), for any rj > 0 satisfying rj < min(/j,, 1 — p), there exists a 
eonstant (depending on [3, but not on t) Cri > 0 such that 


W(t, b) < 


(1 - b)>^-^(e-^ + 


(1 - + b')-()^+'n)db' 

Proof. We first notice that there always exists a function > 0, compactly supported, such 
that 

/3(a) < +gr,{a). 


1 a 


Thereby, for all b' > b and all r > 0, we have 

pe^b' 


B{e^b) - B{e^b') = - ^ /3(s)ds > {n + v)ln - ||gr,||i, 


and 


gB(e^&)-B(e^6') > g-||s,l|i ^g—^ (g”"^ -(- 5)'^+’? 

Then we recall that by definition 

3--B(e"b(l _ 5)M-1 (1 _ 


W{T,b) = ^ 


< 


/g '>')(! - b'y-^db' /j, eB{^-b)-B(e-b')i^i _ yy-idb' 


Inserting (10) in the latter, we obtain the result with C,, = 

It is worth noticing that we can establish with the same proof 

Ve<l/2, / IT < C; 

Jo 

We denote Wref for /3(a) = and notice that in our general case 


□ 


( 10 ) 


□ 


/3(a) = 


X d 


■gia)- 


We denote G{a) = f^g and we can write 

IT(t,6) = K(T)e^^^^'>-^^^^^^Wref(T,b). 

for some K(t) > 0 that insures the normalisation Jq W = 1. We introduce some ?/ > 0 as in 
lemma 1^ We already establish in the proof of lemma [l0| that for £ < 1/2, 

W < 

Therefore, 

1 > K(t) b)db > 1 - 

Let £ > 0 be fixed. Furthermore, since g{a) = o(^), we have 

sup |G'(e'^) — G'(e'^£)| = o [ [ — ] = o (ln£) = o(l). 

b>e yje'^s ® / 
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In particular, for e > 0 fixed, 

gG(e")-G(e"&) ^ 

uniformly on (e, 1) as r —)■ +c». As a consequence, we have for all e > 0 


^ eG(e")-G(e"b)^^^^(^^5)rf5_^ Wref{T,b)db-^0. 


This leads, for any e > 0, to the bounds: 


1 1 - C 

limsupATfr) < -, liminf ATlr) < - j -, 

+- L'wJ +~ J 


Letting £ —>■ 0, we obtain 
What we established proves that, for any e > 0 


lim iL(r) = 1. 
+ 00 


We can conclude using lemma that we have 

[ |W-W,e/|^0. 
JO 


Consequently, 


lim||W-Woolli =0. 
+00 


This ends the proof of lemma □ 

The main property of the pseudo equilibrium is the following 
Proposition 9. W satisfies the following system : 

j drW{T, b) + 5h((l - 6)W(r, b)) + e^/3(e^6)W(T, b) = W(r, b)C{T)S{T) 

{ W{t, 0) (1 + S{t)) = /o' e-/3{e^b)W{T, b)db, 
where S(t) is defined by the equation 

^ C'jr) _ 

^ > {C{r)y C{t)- 

Proof. By computing the partial derivatives of W with respect to b and to r, we obtain: 

96((1 - b)W{T, b)) = -W{t, b) [m + (1 - &)e"/3(e"5)] 

and: 


□ 


( 11 ) 


( 12 ) 


drW (r, b) = W(t, b) 


C'jr) 

C(t) 


— be^l3{e'^b) 


Therefore, W satisfies : 

drW{T, b) + db{{l - b)W{T, b)) + e^/3{e^6)W(T, b) = W{t, b) 


C'jr) 

C{r) 




If we take into account that Vr > 0 ||IT(r,/H^i = 1, by integrating the previous equation 
over b G [0,1], we obtain the value of W(t, 0), hence the claimed system. □ □ 
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The next results justify that 0 is close to Q. 
Lemma 10. Under hypothesis (HI), we have 


lim C{t)S{t) = 0 
•—>■ + 00 


where 


C(t)S(t) = - A* = / [6e'"/3(e^6) - /i] W{t, b)db. 

b'[T) Jo 

Proof. We recall first, by definition 


( 13 ) 


|C'(r)5(r)| < / \e^bfi{e'^b) — p,\W{T,b)db. 

Jo 

We can split then the integral into two parts 

|C'(r)(5(T)| < / \e^bl3{e^b) — fj,\W{T,b)db+ / \e^bf3{e'^b) — p,\W{T,b)db. 

Jo Je-^/2 


b(T) 




Firstly, we have 

To estimate Ii{t) we notice 


I 2 {t) < sup |a/3(a) — ^1 —>■ 0, r —>■+oo. 

(e'^/2 ,+oo) 


p —/2 


/i(r) < sup |o/3(a) - ^1 / IF < sup |a/3(a) - ^| ( ||W - Woo||i + / Wc 
IR+ Jo K+ y Jo 

We already know from lemma II ~ bFoo||i — >■ 0. Furthermore, for large t we have 

c e“'^F-M)/2 


R —/2 


0 < 


VFoo < 


^ 0 . 


Jo (1 ~ 1 — M 

Therefore, we have lim+oo L = 0, which concludes the proof of the lemma. 


□ 


□ 


2.3 Dissipation of entropy with respect to W 

We now introduce the most important tool we will use: the relative entropy (similar to the 
entropy rate of a stochastic process, or the general relative entropy used in I16j i. 

Definition 11. Let w be a solution of the equation Q with support included in [0,1). Let H 
be a convex, continuous function, by parts, which reaches its minimum, 0, at 1. We define 
the generalised relative entropy as: 

And for a non-negative measure n on [0,1) the entropy dissipation DJI(uliy) is defined by 
DH(uliy) = H{u{b))dv{b) u{b)dv{b'^ . 
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Note that DH{u\v) > 0 if is a probability (by Jensen’s inequality). 

We are now in position to establish a first important inequality on the relative entropy 

Proposition 12. Under (HI), the entropy % satisfies the following equality: 

n'ir) = -C{T)DH{u\djr) + C(t)S(t) [ {H{u) - uH\u)) W{t, b)db (15) 

Jo 

where d'jrib) = ^ ® non-negative measure of mass (1 + (5(t)) and u = w/W. 

Proof. The mass of dj is immediately derived from the equation on W(t, 0). 

We have then: 

J drW + i9b((l — b)w) + e^fi{e^b)w = 0 

[Srfb + i9&((l — b)W) + l3{e^b)W = W{t, b)C{T)S{T). 

Denoting u = w/W we arrive at 

drU + (1 — b)dbU = —CSu. 


We multiply this equation by H'{u) and get 


dr{H{u)) + (1 - h)db{H{u)) = -C5uH'{u). 

IWdt(ff(u)) + (1 - b)Wdb(ff(u)) = -WCSuH'(u) 

[HdrW + (1 - h)HdbW + (e^/?(e^6) - 1) HW = C5HW 
dr{H{u)W) + i96((1 — b)H{u)W) + j3{e^b)H{u)W = WC5 [H{u) — uH'{u)]. 

'-V-' 

V 

Taking the integral over b of the previous expression yields: 

^ / H{u)Wdb — W{t,0)H{u{t,0)) /3{e^b)H{u)Wdb = [ rjdb 
Jo Jo Jo 

and finally, 


A J‘ . dr) [h (/ - J' mu) 


yi3{e^b)W 

W{t, 0) 


db 


rjdb 


which, by definition of the entropy dissipation, proves the proposition. 


□ 


□ 


Remark. C{t)5{t) appears naturally as a remainder we will have to estimate in order to prove 
convergence-related properties for Pi, and also for \\W — Woo||i. 


2.4 convergence (withont a rate) to W 

In this section we prove proposition]^ We take therefore H{x) = \x — 1\, we have then. 


DH{u\d'yr) 



\u - l\d"fr - 



ud'Jr — 1 


The core of the proof is the following 
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Lemma 13. Under hypothesis (HI), we have 


\w-W\db<-e^f3{e^) [ \w - W\db + 2\C{t)6{t)\. 
dT Jo Jo 

Proof. By definition, 

l3{e'^b)W{T,b) e'^/3(e'^), 


( 16 ) 


d'Jr = 


W{t,Q) - C{t) 
Furthermore, since /^(u — 1)W = w — W = 0, we have 

DJI(ujdjT) = f |m — 1M7t — f ud'-fr — 1 

Jo Jo 

= f |m— 1 |c? 7 t— f {u — l)djT — S 

Jo Jo 

= [ |m—11^7,-— f {u — 1) {d^r — K{t)W)—6 
Jo Jo '' V ^ 


W{T,b) = K{T)W{T,b). 


>0 


= K{t) f |m —1|1F+ / \u — \\{d'yr — K{t)W) — f {u — l){d'yT — K{t)W) 
Jo Jo Jo 

>K{t) f |m — 1|VF — |i5(t)|. 

Jo 


-|(5| 


Since we also have, for Fl(a;) = |x — 1|, 

\H{u) — uH'{u)\ = I — sign(u — 1)| < 1, 


we obtain for this case 


{H{u) - uH'{u))W 


< 1 . 


And since C{t)K{t) = e'^/?(e'^), we obtain equation (16|. 


□ 


□ 


Since we already have by hypothesis e^l3{e^) —>■ /r > 0, standard ODE arguments yield 


lim sup 

+CX0 


\w — W\db < 


21imsup+^ \C{t)6\ 


We can conclude the proof of proposition using lemma 

Remark. We let the reader check that, by defining /? as in (§, we may replace the non-increasing 
/3 hypothesis by aj3{a) - p, obtaining the following equation instead of (16): 


— a—>-oo 


d_ 

dr 



W\db < -e^^(e^) [ 
Jo 


W|d6-b2|C'(r)(5(T)|. 


(17) 


Remark. We have now finished developing a framework which allows us to deduce the behaviour 
of the entropy R from suitable hypotheses made on /3, and showed that, under mild conditions, 
the entropy tends to 0. The following section will extract a convergence rate from more restrictive 
hypotheses. 
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3 Rates of convergence to the pseudo equilibrium 
3.1 The key situation: /3(a) = 

We consider it best to start by presenting this simple case, since the following proofs contain the 
key innovative elements of the general case while simplifying the presentation of our results. 
Here, the pseudo equilibrium W becomes, with c(r) = e~^'^C(T): 


W(t, b) = 


C(t) 


c(t) 


(1 - 6 ) 1-^(1 + (1 - b)^-f^(e-^ + b)f^' 

We can now compute S as follows. 

Lemma 14. 


S(t) = - 


1 + e- 


(18) 


Proof. 


= Ibe^^(e^b) - /r] W{t, b)db = ^ [5 - (e"" + b)] W(r, b)db 

= —e~'^J^e'^(3{e^b)W{T,b)db. 

By applying Proposition!^ we obtain: 

C{t)S{t) = -e-^(l + (5(T))W(T,0) 

= -e-^(l + (5(T))C'(r) 

resulting in the claimed equality since C doesn’t vanish. □ □ 

Definition 15. We call Coo the limit at oo of c(t) when such limit exists. Here, it is easy to see 
c is a decreasing function and Coo is well defined. 


Lemma 16. 


_ < C'(t)<5(t) < < 0. (19) 


Proof. We have c(r) = e thus establish two trivial bounds on 

c(t). By definition: 


d(r) 


1 


1 


c(t) Jq (e“’‘ + 6)^(1 — 6)1“^ 
hence, since r > 0 and 0 < p, < 1, we have: 


db 


> /(I +6)-»(i - ir-M. > |(^)' d. > 1 1(1 - .,.d. = ^ 

/•i 1 

< / &-'^(1-6)^-M6 = — . 

Jo 


1 

n) - 4 


and likewise: 


1 

c(t) 


It follows that: 

0 < Coo < c(r) < 2(1 + n) < d. 
The result of Lemma [14] allows us to conclude. 


□ 


( 20 ) 

□ 
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Putting together this lemma and equation (16) gives us 


Corollary 17. The following inequality holds 

n'{T) H{r) + ( 21 ) 

All is ready to prove Theorem by applying Gronwall’s Lemma to the previous inequality. 

Proof of Theorem\^ We set f{s) = — which gives fj f(s)ds = —piln ^ Corollary 

[T7| implies that: 


dr 


exp (m +/(s))ds)'H(r) < 8 e exp (^ +/(s))d; 


By integrating over r, we obtain: 

■H(t) < 'HiO)e-fo i^+f + 8 /^+fe-fo z^+ldr' 

< e-i^^n(0)e- fo f + 8e-^^ fj ^dr' 

< e-'^"H(0) (i^)'' + 8e-'^"/;e(2M-iW (i^)'‘dr' 
(T+i=r)''[H(0) + 8/;e(2A^-i)"'dr'' . 

□ 




□ 


3.2 A larger class of /3 


We now consider /3 satisfying (H2). 

Remark. The bound on G is not necessary to prove lemmas 18 and 19 5 € is a strong enough 
hypothesis. However, that precise bound in necessary for our convergence rate estimates and as 
such, we assume it holds throughout the section. 

We take the following notations: 


G{a) 

Wrefir, b) 


W{t, b) 


ir l5'(s)|ds 

Cr-.fjr) 

C(t) exp(- g(s)ds) 

(l-h)i-('(l+e^ 6 )A‘ 


where Gref and C ensure ||lPr.e/(T)|j 1 = ||lP(r)|| = 1 . (We use the same notation as in the proof 
of lemma . 

We have: 


G{t)6{t) 


J^[e^bp{e^b)-fi]WiT,b)db 


Jo e'^bg(e'^b)Wdb + /r fg 


e^b 

1 +e^h 


1 1 T(t, 6 )d 6 . 


Lemma 18. Assume (H2) holds. Then there exists K > 0 such that for any r > 0 and 


0<b<l: 


(1 +e^ 6 )'"(l - 6 ) 1 -/" 


< W (r, b) < 


ATe'"^ 

(1 +e^ 6 )^(l - by-t^' 


( 22 ) 
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Proof. We have: 


W{t, b) C{t) 


Wref{T,b) Cre fir) 


exp 




Cref{T) 


And since W(r, b)db = Wref(T, b)db = 1, it follows that: 


.-llslli < < ell 


Crefi-r) 


which gives us G L°° in the sense given above. 


□ 


□ 


This result leads, through a proof analogous to that of lemma to the following 
Lemma 19. Under hypothesis (H2), there exists a positive M such that: 




e^b 


- 1 


W(r, b)db 




(23) 


_1 + e^6 

We now give the strategy for estimating the rate of convergence. It is based on the same pro¬ 
cedure as before. Consider a non-increasing /? and equation (16), which measures the dissipation 
of entropy (or equation 0 under the corresponding hypothesis). 

n'{T) < -e"^(e")H(T) + 2 |C(t)<5(t)|. 

We have then 

(t) < 2|C(r)(5(r)|e®(""). (24) 


We recall 


C'(t)(5(t) = f {e^b/3{e^b) — p)W{T,b)db. 

Jo 


And also by the definition of W and lemma 


K 


W(T,b) < 7- 

^ ^ “ (e-^ + by{l - 6)1-^ 


Therefore, we easily obtain 


\CiT)6{T)\<K 


e'^b\g\{e^b) 


-db- 


Lemma 19 then gives us 


0 + b)f^{l-by-t^ 

^ e'^b\g\{e^b) 


U 


|C(r)<5(r)| < K 


e^b 

iPe^b 


db + Me^^-^'^^] . 


- 1 W(r, b)db 


)■ 


/o (e-^ + bY{l-bY->^ 

The constant K may change value from line to line. 

We integrate equation (24) and get 

<'HfO')-I- a: f [ 

n(r)e - by 

Using the fact the B{e^) — pLT is bounded from above and below, we can replace B{e^) by pr 
with just a change of constants. 


e^ b\g\{e^ b) 


7— dbdr' + Kf 

Jo 


n{T)e>^^ < K ll + y 


e^^ e^ b\g\{e^ b) 
(e-^' + b)i^{l - by-f^ 


dbdr' + J dr' j . (25) 
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Remark. We let the reader check that the non-increasing /3 hypothesis may be replaced by the 
following condition on the function defined in (Ibl): o/3(a)- > la. Replacing the use of B by that 

■ — n —Vnn 


of ii(a) = Jq /3, we still obtain equation (25) up to multiplication by a constant, since B_{e^) — /xr 
is also bounded from above and below. 

Now the work is focused on the estimate of the middle quantity 

CT rl QtJ-r qT b) 




- dbctr' 


10 JO (e-"'+6)^(1 

We integrate by parts with respect to r' and recall that -^Gie^b) = —e^h\g{e^h)\. We have 


I{t) = 


-^(e^ b )7 -;--TTT- 


db 




lo lo \ [e-^'+by{i-by-t^ ^ {e-^'+by+'^{i-by-t^^ 

The first term is bounded from above (since G > 0) by 

/■i ii„ii 


dr'db. 


Finally, since < 1, 


/o [i + byiy-by-i^' 


1 />r 


-db 


I{T)<K + 2g / G{e^ b). -;-——-— . 

io ^(e-'+6)^(1-6)1-/^ 

We need a sharp estimate on the second term. We focus our efforts on the case 

K 


G{a) < 


[l + aY 


for some a > 0. 


Proposition 20. Assume hypothesis (H2). 
Then, if p, Y 1/2 


nr) < 


K{{l + T)e-^^y, 
R:(e-“^ -he-^^). 


If g = 1/2, then 


nr) < 


K{e--n, 

K{{1 T^)e“’'/^) 


if pL = a<l — pL 
if p,ya<l — fj, 
if a > 1 — pt 
if a = 1 — pL. 

ifa< 1/2 
if a > 1/2 
if a = 1/2. 
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Note that for a > {1 — fj.) the rate is the same than the one for <7 = 0. 


Proof. We need an estimate of 

i-l I-T 2^ rl rT gC/i-a)!-' 

Jo Jo (1 + (e-^' + 6)^(1 - 6)1-^'^^ “ X X 


Putting all together with (25), we have 

/"I P'T 

H{t) < K ( e-'^^+e"'^^ 


al/i-a)!-' 


/o ^0 (e-'+ 6 )A‘+“(l- 6 ) 

We essentially need to estimate the middle term. 

Lemma 21. Under (H2), the following holds true: 

' K fJ e('"-“)^'dr' , 

1 rr g(M-a)-r' 


—dr'db + e-f^^J^ . 


0 Jo ^(e-^'+6)^+“(l-6)i-^ 


dr'db < < 


ifa<l-n, 
K fJ dr' , if a > I — fj,, 

K fj(l + dr' , if a = 1 — /j,. 


Proof of the lemma. Case 1: a < 1 — /i. 
For a < 1 — /r, we simply use the fact that 


1 


db < 


-db < +00 


Jq (e-^'+ 5)'"+“(1 - 6 ) 1 -'^ ~ Jq 6 m+«( 1 -&)i-m 

Case 2: a > 1 — fi. 

It is obvious that we can restrict to the case a < 1. We first need a few intermediate computations. 
Firstly, for 7 > 0, we have 


7 


-- [ —‘^f{b)e^^^^db. 

1 + e-^io 


Where 


/o (e-^+ 6)'y+i(l - 6)1-^ 

f{b) = — 7 ln(e“'^ + 5) + 7 ln(l — b). 


Therefore, we have 

Jq (e^ + 6 )')'+i(l - 6 ) 1-7 “ 1 + e- 


7 


1 


e-^ + 6 


1 + e- 


(26) 


Using this computation and noticing 7 = /i + a — 1, we can easily establish that, for a < 1, 


1 


-db = 


(1-6) 


1 — a 


Jo (e-^ + 6)1‘+“(1 - 6)1-'" 

Applying then (26) and using (1 — 6)1““ < 1, we arrive at 

/■' 1 


(e-^ + 6)T'+i(l-6)1-7 


db. 


lo (e-^ + 6 )'"+“(l - 6 ) 1 -^" 


db < Ke^\ 
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Injecting, we obtain 


Ia{T) <k[ dr' = K [ 

Jo Jo 


Case 3: a = 1 — ^. 
In this case 


rr ,1 g(2M-l)T' 


(2m-1)t' 


+ 6)(1 - by-t^ 

Cutting the integral on & at 1/2 for instance, it is easy to establish 

CT pi 




thereby, we have 




I 


(e—r' + h) 


dbdr^ < K (^j ^I + log(I + e'^ dr'^ 


/a(r) <K f (2 + dr'. 

Jo 


This ends the proof of lemma 


□ 


□ 


To end the proof of proposition!^ we essentially just need to discuss whether the integrals 
of type fj dr' take value or r, and similarly for integrals of type /J" r'e^'^ dr'. □ □ 

4 Rates of convergence towards the equilibrium W^oo. 

4.1 Quantitative estimate of ||W^ — W/x)||i 

In what follows, we justify how the rate of convergence of w to IT can be extended to quantify 
(up to a multiplicative constant) the rate of convergence towards IToo. The main remark is the 
following. 

Lemma 22. Under hypothesis (H2), we have 

llW-IToolli <2 [ le^bfdie^b) - p\W{T,b)db. 

Jo 

Proof. We already know from lemmalim+oo ||lC('r, •) — Woo||i = 0. Therefore 

liW^(T,-)-il"oo||i < J |IT(T',6)-IToo(6)|d6 

Since we have 


/ OO 

J \drW {t , b)\dbdT'. 


9rW = - e"6/3(e"6)IT = -p]W+{f,- e"6/3(e"&))IT, 


it follows that 


And since 


we can conclude. 


\drW\ < 


C\r) 


C'{r) 


C{r) 


M 


+ [ \e^b(3{e^b)-yw. 
Jo 


C{t) 


< [ \e^bl3{e^b) - p\w, 

Jo 


□ 


□ 


18 















Remark. The bound on G has not been used in the proof of the previous lemma, for which g G 
is a strong enough hypothesis. 


We encounter yet again the quantity |e'^ 6 /?(e'^&) — g\W, for which we have already given 
a time-weighted average estimate in the reasoning following equation (25). Let us now provide 
a pointwise estimate. In the reference case, we have 


roo 

\\W{t, ■)-W^\\i<K — - 

Jr 1 + e ^ 

In the situation described by proposition 20 with /3 = -|- g{a) and for some a > 0, we have, 



K 

(I 


In this case, we can split |e'^6/3(e'^6) — /r|IT into two parts and use the previous arguments of 
the proof of lemma [T^ to claim 


roo nl 

/ le^bPie^b) - g\W < / / 

Jo Jr Jo 


gh 




1 


+K 


KMe" 6)17-^ 


(e-^' -f 6)^(1 - 6)1-'" 

1 


(e-^' -f 6)'"(1 - 6 ) 1 -'" 


dbdr' 

dbdr' 


The first term is already know to be bounded by i^"" by lemma 16 . The second term 

satisfies 





< 


< 


< 

< 


/*oo /*! 


/r Jo 
pi 


|e^ bg{e^ 6)| 

1 


1 


6'‘(1 - 6)1 '" V7e-b 


6'‘(1 - 6 ) 1 -'" 
nOO \ 

Ifll I 


d 6 dr' 


1 


1 


Jo 6'‘(1 — 6)1-'" (1-I-e"'6)' 


-db 


Hence the rate of convergence ||IH(r, •) — IHoo||i < K (e^'" i^"^ -|- e 


4.2 Possible influence of the initial condition 

Let us prove a lower bound on the convergence rate of In (|| w(r, ■) — W (t, Ollibo.i)) for initial 
age distribution w{0, 6 ) = 60 ( 6 ). 

Proposition 23. Consider the reference case 

- TT^' 

Suppose the initial age distribution satisfies: 

w°{b) = 6o{b). 

We can bound below the total variation: 

||IT(r,6)-u.(T,6)||Tv >6-1*". (27) 
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Proof. For 


0 (t) 


we have: 


dr 




It follows that 


= 1 - e-^ 

= exp(/J^e'*/3(e'*(;!)(s))ds) 




+ drW + (p'dhW 

[OrW + db ((1 — 0)w) + e'^/3{e’’(/))w] = 0. 


w{t, (fir)) = exp J I3{e^(j){s))ds^ w{0, 4>{0)) 


which, after injecting the corresponding values, yields 

ic(t, 1 — e”"^) = e~^'^S 


with S a Dirac mass. Therefore : 


w{t, ■) = e + Wj 

where Wj > 0 is the distribution of particles that have jumped at least once over (0,r]. 

Since W{T,b)db is an atomless measure, any Dirac mass and Wdb are stranger measures, 
hence: 

\\W{T,b) - w(T,b)\\Tv > 

□ □ 

For fj. < ^, this lower bound agrees up to multiplication by a constant with the upper bound 
given in theoremour convergence exponent is optimal for /i < 

Remark. It is worth noting that we have the trivial bound: 

•H(O) < 2. (28) 

('H(O) can be greater than 1 if w^{b)db has atoms.) 

Remark. Our results are proved for compactly-supported initial age distributions, and they will 
most likely hold for initial age distributions that decrease fast enough. However, if this is not 
the case, the convergence rates might be affected in a way left for future investigation. 


4.3 Convergence rates for natural variables 

We recall: 


where 


n(t, a) = e b) 



ln(l -I- t) 

a 

1 +i ■ 


Definition 24. We set: 


N{t, a) = e b) 


c(ln(l -I- f)) 

(1 -I- t)^(l -I- a)^(l + t — 


which leads to the following Proposition. 


( 29 ) 
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Proposition 25. If 


||t.(T, •) - WiT, OIIlhio.U) < 

then: 

ll"«, ■) - iV(i, (30) 

Proof. The e~'^ appearing as the Jacobian of the change of integration variables is compensated 
by the e~'^ in the definition of w and we get the claimed result. □ □ 

Therefore, in the reference case /3(a) = the distribution of walkers that have age a at 
time t converges to N(t, a) algebraically fast, with a rate that is essentially given by 
we recover corollary 


5 Monte Carlo simulations 

In order to illustrate the evolution of the age distribution of the system and check the accuracy 
of the convergence rates to self-similar equilibrium, we have carried out Monte-Carlo simulations 
for our reference case /3(a) = 

In these simulations, we describe explicitly each individual walker i by associating it with an 
age Qi and a first jumping time r^. The initial age of each walker is chosen according to some initial 
distribution, for instance uniform distribution in [0,1] or a Dirac distribution at age a = 0. The 
first jumping time of each random walker is sampled from the distribution ^(r) = 
that corresponds to our reference jump rate /3(a) = p,/(I -I- a). The simulation then iterates 
the following steps: (i) find k, the walker with the earlier jump time: k = argminr^, then (ii) 

i 

make it jump, i.e. reset its age a^ = 0 and finally, (Hi) pick its next jump time according 
to ^(t). During the simulation, we store the distance between the dynamic equilibrium W at 
that time and the observed distribution of rescaled ages bi = ai/{l + t) of all the walkers i in 
the simulation: ||r(;(r, •) — kb(r,We also compute at each time step the norm 
of the difference w{t, •) — Woo- Unless stated otherwise we use 20,000 random walkers in each 
simulation. 

First, we note that in all cases, the simulated distance between w and the pseudo¬ 
equilibrium W is indeed bounded above by the expression given in theorem (except at very 
high T, when our bound becomes lower than the numerical error of the simulation). The example 
given in figure corresponds to /i = 0,4 and p, = 0,8, and an initial age distribution Dirac at 0 
for the red dots, and uniform on [0,1) for the blue dots, the black curve representing the upper 
bound proved in theoremtaken for 33(0) = 2, which is an upper bound for 33(0). As we see, 
the multiplicative constant we lose (the overestimation of K in theorem corresponding to the 
losses throughout the inequalities used to prove our bound) is not too high. 
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Figure 1: ||u;(r, •) — VF(r, •)||i lies under the theoretical bound (black curve) for an initial 
age distribution 6o (red dots) and Z^(0,1) (blue dots). 


In order to illustrate graphically the behaviour of the solution to our equations and its 
convergence towards the pseudo-equilibrium, Figure displays, for /i = 0,6 and an initial age 
distribution = Sg, the time evolution of the simulation results expressed either in 

the original variables n{t,-) (histograms), (full line) on the left-hand side column or in 

the rescaled variables w{t,-) (histograms), IF(t, •) (full line) on the right-hand side column. 
Moreover, the rescaled variables panels also show as grey dotted lines the equilibrium Woo, to 
which W converges. 
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Figure 2: Evolution of n, N, w, W and Woo along time, for = 0, 6 and an initial 
distribution nP = = 6o. 


age 


From visual inspection of the is figure, it is clear that n(t, a) largely flattens as t —> oo (note 
the difference in the y-axis scale between the panels). The figure depicts a pointwise convergence 
of the simulated w to the pseudo-equilibrium W which in turn converges pointwise to Woo- 
Moreover, it illustrates how rescaling allows a better description of the self-similar behaviour, 
which is difficult to grasp in natural variables since n converges pointwise to 0. The next sections 
quantify the simulated convergence rates. 
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5.1 Exponential fit of ||ty(r, •) — lE(r, ■)||i 

To quantify the convergence rates in the simulations, we fit the distance ||w(t, •) — W{t, •)|| by 
the following function: 


/(A,A,B) =y4e-^^ + Be-(i-^)^ + C'. (31) 

Remark (Heuristic estimate of the error term). C is a simulation error, that we evaluate to C « 
0,1. This is consistent both with empirical evidence and with a simple heuristic overevaluation 
of C as =ffhins/tides, which is roughly 0,16. 

Remark. According to the above analysis one expects X = fi. A and B are multiplicative 
parameters: we expect A around 2 and |i?| close to 0, since 'H(O) = 2 and our upper boundary 
is of the form 'H(r) < ['H(O) — (< 0)] e~^'^ + ke~^^~^'''^, with k small. 

Remark. Another possible explanation of the predominance of e~^'^ over in the conver¬ 

gence rate is linked to the fact that, for a given /3, two solutions wi and W 2 corresponding to 
different, compactly supported initial conditions, satisfy, for a certain constant K (see Subsection 

lk(r,-) - u>i(r,-)||i < Ke 

Figure [^presents as examples three cases that exhibit a certain diversity : ^ = 0, 9, // = 0, 5 
and /r = 0,2. We plot in red dots the evolution along r of the simulated value of ||w(r, •) — 
1F(t, •)||i and use function / defined in equation to fit the results (blue curves). The fit 
results are given in the figures, ± one standard deviation. 

We first note that in the three panels of figure]^ C « 0.1 as expected and our estimates for 
A are very close to p-. Note that in the second panel, with p = 0.5, the values of A and B cannot 
be estimated independently thus the large inaccuracy/variance on their determination. Finally, 
the third panel shows a marked discontinuity around t = 6. This is due to the discretisation 
of the age distribution: with small values of p, the number of random walkers that have never 
experienced a single renewal during the simulation period becomes large. Since, according to our 
initial conditions, all walkers have the same initial age, many walkers will enter the last age bin 
simultaneously thus causing the observed discontinuity. However even in this case, we obtain a 
very good fit for A by restricting the fit to the values before the discontinuity and fixing C to 
0 . 8 . 

Figure 1^ summarizes the values of A determined from Monte-Carlo simulations identical to 
those shown in Figj^ (red crosses), together with the diagonal line A = p (blue). For all the 
values of p tested, the simulations confirm that w tends to W with a sum of exponential rates 
given by p and 1 —p. Therefore, taken together, those simulation results, while agreeing with our 
analytical estimations, suggest that our estimate of Ijw — VF||i may not be optimal, in particular 
for larger values of p. 
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A = 0.899 ± 0.007 
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A = 0.523 ±0.033 
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Figure 3: Fit by / defined in equation ( [^ (blue curves), for different ji, of the simulated 
||t(;(r, •) — VF(r, •)111 (red dots). Initial age distribution: = 5q. 












Figure 4: Values of the exponent A found by using function / from equation (31) to fit 
the simulated values of \\w{t, •) — W{t, •)||i for /r G {0,1; 0, 2 ; ... ; 0, 9}. 


5.2 For large /i, W provides a better asymptotic approximation of w 
than Woo 

Figurecompares the distances between w and W (red dots), w and Woo (black dots), or W 
and Woo (green curve), for three values of /r. This figure shows that for for ^ < 0.5, W{t, •) and 
Woo are systematically much closer to each other than to w. However, as /i increases, this trend 
reverses: for large enough t, w becomes significantly closer to W than to Woo: the distance 
between w and W converges much faster. Therefore, according to those simulation results W is 
a much better asymptotic approximation w for fi > 0.5, thus justifying further its utility here. 


6 Future developments 

Throughout the article we have estimated norms, but we have presented the estimates in 
the context of an entropic structure. It is indeed possible by means analogous to ours to prove 
entropy inequalities for dissipations corresponding to other H functions than 11 — • |. For instance, 
the classical H(x) = xhix — x + \ also allows us to prove a convergence rate of the corresponding 
entropy to 0: it is also Thanks to the Csiszar-Kullback inequality, it is also 

possible to prove a rate of convergence of ||w(t, •) — W(r, •)|ji to 0, albeit one worse than that 
obtained in theorems |4] and [5l 

We may encounter inequalities such as that of proposition bounding the derivative of 
an entropy with respect to a probability measure Wdb by an entropy dissipation with respect 
to another measure (which we can compare to the dissipation with respect to a probability 
measure dy,-). When the comparison of DH{u\Wdh) and DH{u\d')) does not follow calculations 
as straightforward as ours, an alternative may be to rely on a precise Jensen estimate comparing 
the entropy dissipations with respect to two absolutely continuous probability measures. 

Here, we have considered a spatially-homogeneous (zero-dimensional), age-dependent renewal 
probability /3(a). We believe the ideas we have exposed may be used to tackle the problem with 
a spatial extension, for instance in a discrete space setting. 
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Figure 5: Influence of fj. on ln||r(;(r, •) — IF(t, •)||i (red dots) and ln||t(; — IFooH (black 
dots): for higher values of fj,, w is signihcantly closer to W than to VFoo- 
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One major interest of our age-structure approach of CTRW is that the dynamics remain 
Markovian. We believe that keeping Markovian properties will be crucially helpful when intro¬ 
ducing the coupling between sub-diffusive CTRW and reaction, since the coupling should simply 
consist in the addition of the reaction and the subdiffusion terms (contrarily to the case of frac¬ 
tional dynamics). However, the extent to which the supplementary age variable will make this 
process more complex remains to be evaluated. 


7 Appendix 


The case /t = 1 

It is quite interesting to notice that even if the behaviour is not really self similar, our method 
gives a precise asymptotic for the case ^ = 1. To illustrate this, we focus on the reference case: 
/3(®) = T^- this case the ’pseudo equilibrium reads’ 

W (r, b) = 


(e-^ -f b) log(l -f e^) ■ 

This pseudo equilibrium tends to a Dirac mass but gives still quantitative information. Indeed, 


following the same computation than for equation (16) for the case /r < 1, we obtain easily 

d '■i '■1 

d-T JQ r -r e .70 

Where we also have 


\w-W\< f\u,-W\+2 \C{t)S{t )\. 

Jo -L -r e Jo 




This leads to 


C{t)S{t) = - 


1 


And finally. 


0 (e“^-I-log(l-I-e"’") log(l-I-e"’") \e“'’'-I-1 e 
C{t)S{t) = ———^—7^-—^ —>■ 0. 


(1 -I- e"’’) log(l -I- e'’") 

And we can still claim that \w — W\ —>■ 0. We can give a (rough) estimate for a rate of 
convergence. Integrating, we have 


/o 




We estimate the second term 
1 / 


e'’’ 1 

-T^dr' = 


1 


-du 


l + e^jQ log(l-l-e^') l + e'^ log(l-f u) 

This term behaves as r. Indeed, we have easily (splitting the integral at e"’’’ for a < 1. 


1 1 

< 


du < 


log(l-I-e’’) 1 + e'^Ji log(l-l-u) l-l-e”'’" log(l-I-e"’’) 


f \w-W\< ^ ^ f \w - W\{t = 0) + 

Jo + 6 JO i- + T 


K K' 
< 


+ T l + T 
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